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We investigate a classical lattice system with N particles. The potential energy V of the scalar 
displacements is chosen as a <f> 4 on-site potential plus interactions. Its stationary points are so- 
lutions of a coupled set of nonlinear equations. Starting with Aubry's anti-continuum limit it is 
easy to establish a one-to-one correspondence between the stationary points of V and symbolic 
sequences a — (<ti, . . . ,ctjv) with a„ — +,0, — . We prove that this correspondence remains valid 
for interactions with a coupling constant e below a critical value e c and that it allows the use of a 
"thermodynamic" formalism to calculate statistical properties of the so-called "energy landscape" 
of V. This offers an explanation why topological quantities of V may become singular, like in 
phase transitions. Particularly, we find the saddle index distribution is maximum at a saddle index 
n max _ -jyg ^ QT £ ^ £ ^ Furthermore there exists an interval (v* , v max ) in which the saddle index 
n 3 as function of average energy v is analytical in v and it vanishes at v* , above the ground state 
energy v gs , whereas the average saddle index n s as function of energy v is highly nontrivial. It can 
exhibit a singularity at a critical energy v c and it vanishes at v ga , only. Close to v gs , n 3 (v) exhibits 
power law behavior which even holds for noninteracting particles. 

PACS numbers: 



I. INTRODUCTION 



We consider a classical N- particle system in three dimensions with potential energy V(afi, . . . , £n)- Both the 
Newtonian dynamics and the thermodynamical behavior is obtained from the knowledge of the function V. Therefore 
one may ask: What are the characteristic features of V which are crucial for the dynamics and thermodynamics of the 
N-particle system? From a mathematical point of view the stationary points, in case of non-degeneracy also called 
critical or Morse points, yield important information on V. These points are solutions of the set of coupled, nonlinear 
equations: 

dV 

—zr(xi,...,x N ) = ,n = l,...,N . (1) 

These solutions are denoted by x( Q ) = {x^\ . . . , 2$ ), a — 1, 2, ... , Mn, where Mjq is believed to be exponentially 
in N . Here some comments are in order. First, if the system is homogeneous and isotropic any translation, rotation 
and reflection of {a;^} is a stationary point as well, with the same potential energy. Accordingly there is a continuous 
degeneracy of {x^}. Choosing the variables x n , e.g. with respect to the center of mass and fixing the orientation 
one gets rid of this degeneracy. This is assumed in the following. Second, if V is harmonic, i.e. 

N 3 

V{x, , . . . , x N ) = - J2 E M ^>L (2) 

= 1 =1 

Eq. becomes a set of linear equations: 

N 3 

5^^M«„a4 = J n=l,...,N; i = 1,2,3 (3) 

m— 1 j — 1 

In the generic case where det(M^ m ) ^ 0, there is one solution, x n = 0, only. Consequently, a necessary condition 
for exponentially many solutions is the nonlinearity of Eq. (JIJ. Having found all stationary points they can be 
characterized by their saddle index n s (x^) — N s (x^)/(SN). N s (x^) is the number of unstable directions, i.e. 
the number of negative eigenvalues of the Hessian ((d 2 V/dx l n dx J m )(K < - a >). Therefore n s varies between zero and one. 

The role of n s for the dynamics is obvious. If the trajectory x(f ) in configuration space is mostly close to stationary 
points with n s very close to zero, i.e. x(t) is close to local minima of V{x\, . . . , xjv), then the motion will mainly be 
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due to thermal activation, i.e. hopping-like. In contrast, if x(t) is mostly close to stationary points with n s almost 
equal to one, a diffusive-like dynamics may occur. Recent MD-simulations for liquids have determined the average 
saddle index n s as function of temperature T— and as function of the potential energy v per particle^* 3 -. It has been 
found that n s decreases with decreasing temperature and energy and its extrapolation vanishes at a temperature T* 
and an energy v* , respectively. It is interesting that T* is very close to the mode coupling glass transition temperature 
T c , at which a transition from ergodic to nonergodic dynamics takes place. This finding corresponds to the change 
from diffusive- to hopping-like motion, when decreasing the saddle index from one to zero. 

There is evidence that the saddle index also plays a role for equilibrium phase transitions. For smooth, finite range 
and confining potentials V it has been proven^ that a necessary condition for an equilibrium phase transition is a 
singular change of the topology of the manifold 

M N (v) = {(fi, . . .,x N )\V{x u ...,x N )< v}. (4) 

Particularly, if Mn(v, n s ) is the number of stationary points of A4n(v ) with saddle index n s , the Euler characteristic 

37V 

x {v)=^2(-i) N 'M N {v,N./3N) (5) 

iV„=0 

can be singular at a critical energy v c . However, for nonconfining potentials the unattainability of a purely topological 
criterion for the existence of equilibrium phase transition has been claimed^. Besides these interesting results it 
has been found for a two-dimensional <f> -model from numerical computations^, for the exactly solvable mean-field 
XY modeli and mean-field k-trigonometric model^ that v c coincides with 7V" 1 (I/) (T c ), the canonical average of the 
potential energy per particle at the phase transition point T c . However, for a mean-field (/> 4 -model it has been shown 
that v c f N- l iy)(T c ^.. 

Finally we want to mention an interesting result for the saddle index distribution function P/v(n s ). For a binary 
Lennard- Jones system with up to N — 13 particles there is strong evidence that P/v(n s ) is Gaussian with a maximum 
at n™ aa; w 1/3^. However, the relevance of this result for dynamics and thermodynamics is not yet clear. 

This exposition so far has demonstrated that the potential energy surface (PES) characterized by stationary points 
and their saddle indices is of physical importance. For further applications of the energy landscape description the 
reader may consult the textbook by D. J. Walesii. 

It is the motivation of the present paper to derive a relationship between and symbolic sequences er = 
(<7i, . . . , <7 n) where a n takes values from an "alphabet". This will be explained in the next section. Particularly 
the usefulness of this relationship for the calculation of statistical properties of the PES will be shown. To make 
this relation more explicit we will investigate in the 3. section noninteracting particles in an on-site potential. In 
the 4. section it will be explored how the features of the PES change under switching on an interaction between the 
particles. The final section contains a discussion of the results and some conclusions. 



II. DESCRIPTION BY SYMBOLIC SEQUENCES 

For some nonlinear dynamical systems it has been proven (see e.g. Refs^i^) that there is a one-to-one correspon- 
dence between the orbits and symbolic sequences. Therefore, the brilliant observation by Aubryi4 that Eq. Q for the 
one-dimensional Frcnkcl-Kontorova model is identical to the standard map, already provides a link between stationary 
points and symbolic sequences. This link exists for a certain class of one-dimensional models^. But we think that 
this relationship may be more general. Let us choose a simple potential with two degrees of freedom: 

V(x±, X2) = sinxi sina^ (6) 

Its stationary points form a square lattice with lattice constant n/2. The question arises how to label these stationary 
points. The first guess is to count the points as shown in Fig.^,. 

However, there is a drawback since two neighboring points do not possess "neighboring" labels (see e.g. the points 
with labels 28 and 53). In order that the labelling preserves the local arrangement of the points one has to choose an 
"alphabet", which are the integer numbers a n € Z. Then the "sequences" (171,172) provide a labelling for which the 
local properties are preserved (cf. Figure ^,+b). Similarly, the counting of the stationary points of an arbitrary PES 
by a = 1, 2, ... , Mat, as done in the 1. section, is not compatible with the local properties of these points. Although 
the explicit determination of an "alphabet" A and of the one-to-one correspondence between the stationary points 
of an arbitrary PES and symbolic sequences (ci, . . . , o~n) with a n G A in general is not possible we believe, however, 
that such a relationship may exist for certain potential functions V . This will be proven in the third and fourth 
section for a certain class of functions V(a?i, . . . ,xjv). Of course, for a finite system there will be a finite number of 
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FIG. 1: Labelling of the square lattice by (a) natural numbers and (b) by "sequences" (cri^cri) with integer values of a n - 

stationary points which can be labelled by a = 1,2, 3,.... But similarly to the simple model © this will be not the 
appropriate labelling. 

Thus, let us assume that we have found A and the one-to-one mapping: 

<t = (er 1; . . . , <7]v) >-» (xi(er), . . . ,afjv(<x)) = x(er) (7) 

between the symbolic sequences a and the stationary points x^ a '=x(tr). Their potential energy is given by: 

E(a) = V(x(«r)). (8) 

Note that depending on the specific potential energy V the sequence may contain more than N symbols, e.g. cr = 
(<7i, . . . , <73jv). The Hessian at these stationary points is a function of <x, too. Therefore its number of negative 
eigenvalues, N s , and accordingly the saddle index n s is a function of <x, only: 

n,(<r) = N s (a)/(3N). (9) 

In the remainder of this section we will demonstrate the usefulness of the mapping Q for the calculation of some 
properties of the PES. One of the most interesting quantities is the joint probability density Pjsr(v,n s ) of stationary 
points with potential energy v per particle and saddle index n s . Without making use of relation J7J) it can be 
represented as follows: 



P N (v,n s )=Mj [ n^|det(|^)|(nn<5(^)) (10) 

5{v - A^- 1 F(x))5(n s - iy -1 JV s (x)). 

One of the main technical problems for the evaluation of Pn(v, n s ) from Eq. i|10|) is the occurrence of the modulus of 
the determinant of the Hessian. This modulus prevents the use of an integral representation by Grassmann variables. 
For mean-field spin glass models it has been shown that the modulus can be neglected, at least for low energies^. 
Now, making use of Eqs. 0-©, this probability density can also be represented by: 

P N (v,n s ) = N 2 M^ 1 J2s(Nv-E(cr))5(Nn s -N s {cT)) (11) 

cr 

which is rewritten as follows: 

oo oo 

P N (v,n s ) = ^^- J dX J ^ e *WM».-/<A./0] (12) 
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with the "free energy" per particle 

f{\,n)=iN- 1 \nZ{\, l i-N) (13) 

and the "canonical partition function" 

Z{\ V,N)=J2 expH(A£(<r) + (iN s (<t))}. (14) 

This demonstrates that, e.g. Pjsf(v,n s ) can be obtained from a "canonical ensemble" with probability density: 

p(tr) = l e -«W+t*f.W) . (15) 

A can be interpreted as an inverse "temperature" and \i as a "field" acting as a bias on the number of negative 
eigenvalues of the Hessian. 

In the thermodynamic limit N — > oo the saddle point solutions \*(v,n s ), 
fi*(v,n s ) of 

V = |^(A, fi), n s = |^(A, n) (16) 

yield up to a normalization constant: 

P N (v,n s ) ~ e +"»(«.».) (17) 

with the configurational entropy (per particle) of stationary points with energy per particle v and saddle index n s : 

s(v,n s ) = i{\*(v,n s )v + (i*(v,n„)n a - f(X*(v,n a ),fi*(v,n a ))}. (18) 

s(v,n s ) allows the calculation of the saddle index as function of energy. However, there are two possibilities. First, 
we determine the maximum v of s(v,n s ) for n s fixed. This yields the saddle index n s (v) as function of the average 
energy. Second, keeping the energy v fixed the maximum of s(v, n s ) yields the average saddle index n s (v) as function 
of energy. These two functions, n s and n s are not identical, in general. From 1)18(1 we find: 

X*(v,n s (v)) = 0, fi*(v,n s (v)) = (19) 

where Eq. {TfJJ) has been applied. 

oo 

The saddle index distribution Pn&s) — J dvP^f(v,n s ) can be represented as: 



oo 



P N{ns) = ^MZ / rf/1 ^K-/.W] ( 2(), 



2tt 
with 

/o(/x) = /(0, M ). (21) 

For N — > oo the saddle point solution fJ^(n s ) of 

n s = ^(M) (22) 

leads to 

PjvK) ~ e Nso ^ (23) 

with 

so(™s) = i{^o{n s )n s - f (Ho( n s))}- (24) 
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Its maximum position n™ ax is given by: 

^ (< aX ) = • (25) 

This last paragraph has demonstrated the usefulness of the "thermodynamic" formalism based on a canonical 
ensemble l|15|l in the space of symbolic sequences. This fact also allows to understand why there can be a singular 
topological change, e.g. as a function of v. There is a lower critical dimension di ow such that there is a critical 
"temperature" A^ 1 for d > di ow at which /(A,/i) is singular. This singularity then also occurs in the "Legendre 
transform" s(v,n s ) of /(A, /i). 

In the third and fourth section we will show that this "thermodynamic" formalism applied to a certain class of 
lattice models allows an analytical calculation of PN(v,n s ) from which the saddle index distribution and the energy 
dependence of the saddle index can be derived. 



III. NONINTERACTING PARTICLES 

We consider a crystal with N lattice sites. Let x„ be the scalar displacement of particle n from its lattice site with 
lattice vector R n . Then we describe the potential energy by a 4 -model: 

N 

V(x) = ^V (x n ) + eVi(x), e>0 (26) 

71=1 

with the double-well-like on-site potential: 

V (x) = -hx - -x 2 + jx A , h>0 (27) 

and interaction V\{x). Choosing an appropriate energy and length scale any general quartic on-site potential can be 
put into the form of Eq. f2*7| . 

Inspired by Aubry's anti-continuum limit (originally called anti-integrable limit )ILi2i we start with e = 0, i.e. 
neglecting the interaction. Then Eq. Q becomes: 

4 - x n = h . (28) 

If h < h c = 2/(3v / 3) there are three different real roots x n = x an {h) 1 a n = +,0,—. Therefore the stationary points 
are given by 

xo(er) = (x ai , . . . ,x aN ) , (29) 

i.e. we have found the one-to-one correspondence between stationary points and symbolic sequences with an "alphabet" 
A = +, 0, — . Using Eqs. i|5|). H27fl and (|29|l we hnd for the energy of the stationary points 

N 

er„ + e 2 er„] (30) 

71=1 

with 

e a = v a ,e 1 = - u_) , e 2 = -v + ^{v+ + u_) (31) 

v a = Vq{x„). 

These results are obvious, since the stationary points of V for e = are uniquely determined by the stationary points 
of the local potential Vo(x). The number of extrema of Vq(x) determine the "alphabet" A and their energy Vo(x a ) 
yields the coefficients e v ,v — 0,1,2 of Eo{a). Furthermore the function N s (cr) is also easily determined. Since we 
identify the maximum of Vo(x n ) with a n = 0, its absolute and local minimum, respectively, with a n = + and a n = — , 
it is: 



iV s ( CT )=£(l-^). (32) 
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Having found the "alphabet" A, the mapping Q as well as Eq(ct) and N s (cr) we can calculate Pn(v, n s ) as described 
in the previous section. The calculation of the "partition function" (|14ll is easy. One gets from Eq. H13|) : 

/(A, fi) = (Ae + n) + i ln[l + 2e~ i{ - Xe2 -^ cos Ae x ] . (33) 
The saddle point solutions of Eq i|16|) are easily determined leading to Pn{v,ti s ) from Eq. (|17|l with: 

s(v, n s ) = — [n_(u, rt s ) lnn_(«, n s ) + n+(v, n s ) lnn + (v, n s ) + n s lnn s ] (34) 

where: 



n±(v,n s ) = ±(v- - v + ) 1 [(v T - v) + (v - v T )n s }. 
for asymmetric double well, i.e. v + ^ v-. Quite analogously one finds the saddle point /io(n s ): 

1 — n s 



2n s 



leading to: 



so(n s ) 



,lnn s + (l-n s )ln(l-n s )/2] 



(35) 



(36) 



(37) 



for symmetric and asymmetric double well potentials. 

To determine the relation between the saddle index and energy for symmetric double well, i.e. for h = 0, we have 
to realize that n s and v are no longer independent variables. This can easily be seen from Eq. I|30|) . taking e\ = 
into account, which follows from Eq. (|31|l . Then we arrive at 



E(tr) = Nv+ + (v - v+)N s (a) 
where Eqs. (|3*T|l and were used. This yields immediately: 

v(n s ) = v + + (v - v + )n s , 



(38) 



(39) 



and no distinction between n s (v) and n s (v) exists. Because of the dependence of both variables v and n s on each 
other their joint probability density reduces to Pn(u s ) which is proportional to the probability density of v. From 
Eqs. (|25|l and l|36|) we find for symmetric and asymmetric double wells: 



1/3 



(40) 



which is the maximum of so(n s ) (Eq. I|37|l ) and therefore the maximum of the saddle index distribution Pjv(^s)- fw("s) 
is a Gaussian for \n s — n" lax \ = 0(l/y/~N), in agreement with the numerical result for Lennard- Jones clusters^. 

Let us return to the asymmetric double well. The functions n s (v) and n s (v) can be obtained from the maximum 
of s(v,n s ) for fixed n s and fixed v, respectively. As a result we get 



and 



n s (v) 


(«_ -+ 


- v + ) - 2v 






- («- +v+) 




■n + (v,fi s )- 


VQ-V_ 


n_(«,n, s )- 


v -v + 


v_-v + 


v_ - V+ 


n s 




n s 





(41) 



(42) 



with v-, vq, w+ from Eq. (|31() and n±(v, n s ) from Eq. I|35|l . n s is a linear function of the average energy v. It vanishes 
at v* = \{v- + v+) which is above the ground state energy v gs — v+, since we assumed asymmetric double wells. n s 
becomes one for v = Vq, the height of the unstable extremum of Vq(x). n s (v) is presented in Figure^. Eq. (|42l) is 
an implicit one for n s (v), which can not be solved analytically. The numerical result is shown in Figure I^Jd. There 



are two features to be mentioned. First, n s (v) vanishes at the ground state energy v g 



proven analytically that it exhibits a power law behavior in v close to v gs (cf. the inset of Fig. |2t>): 



only. Second, it can be 



n s ( v ) ~ ( - Vgs ) s ° + Q(( - Vgs ) a ° +1 ) + Q(( V Vgs 



(43) 
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FIG. 2: (a) saddle index n s as function of the average energy v, (b) average saddle index n 3 as function of the energy v. Both 
figures are for noninteracting particles (e = 0) with vo = 0, v- = —1 and v+ = —1.8. The inset in Fig.^ shows the comparison 
of the exact result for n s (v) from Eq. 1421 (solid line) with the leading power law from Eq. 1431 (dashed line). 



for — 1 <C 1 with an exponent: 



vo 



> 1 



(44) 



It is interesting that we find for n s (v) (in contrast to n s (v)) a nonanalytical u-dependence close to v gs , despite the 
neglection of interactions. In the next section we will show that this power law also exists in case of interactions. 

So far we have demonstrated that the "thermodynamic" formalism can be applied to get the saddle index distri- 
bution, n s (v) and n s (v). n s (v) and fi s (v) also follow directly from the saddle points A* and /i* (cf. Eq. I19fl . This has 
not been used, because X*(v,n s ) is a rather lengthy expression. The results l|37() . (|40|l concerning the saddle index 
distribution, and (|41|l for n s (v) can be obtained much easier as follows. The number of stationary points characterized 
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by N a , which is the number of er„ of <r equal to a e { — ,0,+}, is given by 



N + \N \NJ. 



(45) 



which for fixed saddle index n s — Nq/N has its maximum weight for N + = iV_ = ^(N — N s ). Substituting this into 
expression (|45|l we get by use of Stirling's formula: 

P N (n s ) ~ e Nso ^ (46) 

with so(n s ) from Eq. 1(37} . The maximum position n™ ax is also obvious. If N™ ax = N™ ax = N™ ax = N/3, expression 
l|4*5|l takes its maximum value. Therefore n™ ax = N™ ax /N = 1/3. Next we use E{cr) from Eq. ((37]]) to calculate 
v(n s ) — N^ 1 (E(cr)) . Using again that the maximum weight of (|45|) . for fixed n s , follows for N + = AT_, we get: 

(J2<rn)=N+-N- = 0, (^a2> =N-N (47) 

n n 

Substituting this result into Eq. (|3~U|) yields: 

v(n s ) = -(»_ + w+) + [v -\( v -+ v +)]n s (48) 

from which Eq. I|41|) is reproduced. The question arises: Which of these features remain valid in the presence of 
interactions? This will be discussed in the next section. 



IV. INTERACTING PARTICLES 



Now we will study the influence of interactions on the results derived in the previous section. V(x) is given by 
Eq. iffifl) and l|77|) where we will assume that the interaction energy V\{x) is (i) at least twice differentiable and 
(ii) does not grow faster than quartic. To investigate the labelling of all stationary points of V(x) for e > we 
adopt the method used by MacKay and Aubry and MacKay and Sepulchre to prove the existence of breathers^ 
and to investigate multistability in networks 20 , respectively, i.e. the implicit function theorem. In the last section we 
have proven for e = that all stationary points x (cr) = x(e = 0,<x) are uniquely labelled by symbolic sequences 
c, c n = +, 0, — , provided h < h c . The eigenvalues A„ (0, <r) of the Hessian (e = 0) are nonzero for all <x and h < h c 
with: 



sign X u=n (0, a) = - 1 + 2< (49) 



This implies that 



for all er. Let us introduce the functions 



OXfiOXui 

dV 

4> n (e,x) = - — (e,x), (51) 
ox n 

which are at least one times differentiable. Note that we have made the e-dependence of V explicit. Then Eq. (JTJ 
reads: 

cf> n (e,x) = 0, n = l,...,N . (52) 
We have found solutions xo(cr) for e = 0: 

cj> n (0,x o (cr)) = 0, n = l,...,N . (53) 



Because the Jacobian: 
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for all cr and < h < h c , we can apply the implicit function theorem^! which guarantees the existence of a neighbor- 
hood U(cr) of (e = 0, x (cr)) 6 M jv+1 , of an open set W(cr) — [0, e c (f)] C M and of functions x n (e, cr), n = 1, . . . , N 
(at least one times differentiable) such that 

(j) n {e,x(e,cr)) = , (55) 
for all e S VF(cr) and all cr. This leads to a nonvanishing critical value 

e c = min e c (cr) > (56) 

such that x(e,a) are stationary points of V(x) for < e < e c . Therefore the one-to-one correspondence between 
these points and symbolic sequences cr with a n e {+, 0, — } holds for < e < e c and h < h c . 

At e c a bifurcation occurs at which stationary points disappear and new ones may be created. Accordingly the 
"alphabet" and/or the mapping between stationary points and symbolic sequences will change at e c . This bifurcation 
is signalled by the vanishing of at least one of the eigenvalues A„(e c ,cr), for one or more sequences cr. Consequently, 
it is 

signA„ (e, cr) = signA^ (0, cr) (57) 

for < e < e c , since a change of sign occurs at e c , only. This has the strong implication that the saddle index of x(e, cr) 
is identical to that of a;(0, cr), which is given by the number TVo of <r n = 0. This result can also be put into other 
words: Although the relative heights of the stationary points will change with e, the topology will remain unchanged. 
By this we mean that a stationary point with N s unstable directions remains a stationary point with N s unstable 
directions for < e < e c . We note, that this result requires condition (ii) for the function V\. In case that V\(x) 
grows faster than quartic, additional stationary points may occur already for arbitrary small values of e. Validity of 
(ii) guarantees that x(e, cr), which are continuously connected to a:(0, cr), are the only solutions of Eq. Q. 

Some general conclusions can be drawn from this result. Since the one-to-one correspondence between the stationary 
points and cr is preserved for < e < e c , the saddle index distribution P/v(n s ) and accordingly n" lax = 1/3 remains 
the same. Although the calculation of n s (v) and n s (v) needs the knowledge of the energy E(cr) of the stationary 
points, one can prove some general properties of both functions without using the explicit form of E(cr). Let us begin 
with n s (v) which is obtained from v(n s ). As already shown in the last section v(n s ) is directly obtained from 

v(n s ) = (NM N (n s ))- 1 ^ E{&) (58) 

N s (a-)=n s N 

where M/v(n s ) is the number of cr with N s (cr) = n s N. n s — yields all local minima of V(x), including the ground 
state. Most of these minima have an energy above the ground state energy Nv gs , i.e. it is: 

v(0)>(NM N )- 1 Nv gs M N (n s )=v gs . (59) 

For n s = 1, i.e. N = N s , there is only one stationary point, the maximum of V(x), with energy E(0, . . . , 0) such that 

v(l) = N- 1 E(0,...,0) = v max (60) 

From this we find that n s (v) vanishes at v* > v gs and becomes one at V = v max . Whether or not n s {v) is always 
monotonically increasing with increasing v is not clear. In contrast to n s (v), the average saddle index n s (v) is nonzero 
for all v > v gs . Without presenting a rigorous prove, let us explain why this should be true. The ground state 
belongs to n s — and is characterized by cr with ct^ s = + or — . This also holds, if it is degenerate. Now, let us 
choose K particles j\, . . . ,jk for which we change cr?* g {+, — } into aj k = 0. The corresponding stationary point 
has n s = K/N. If K = 1, we generate a "defect" with excitation energy €j 1 . Then K "defects" have an energy 

K 

Ejx—jK = Nv gs + e ifc+ (•••)> where (...) is the interaction energy of the defects. It is obvious that Ej x ...j K /N — * v gs 
fc-i 

for K — > oo, N — > oo with K/N — ► 0. "Defect" configurations with fixed energy v = Ej 1 ,,,j K /N will have fluctuating 
K with < K/N. If we choose v — v gs arbitrary small but nonzero the average saddle index n s — K/N will be small, 
but finite, too. Consequently n s can only vanish at v gs . Since there is one stationary point (maximum) with energy 
v m ax, only, it must be n s (v max ) = 1, i.e. 

n s (v = v max ) = n s (v max ) = 1. (61) 
Again, it is not obvious whether n s (w) is monotonously increasing with increasing v. 
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To get more quantitative results for n s (v) and n s (v), we have to specify E(a). Since the solutions x(e, a) are 
not known exactly this can only be done approximately. In the following we will present the crucial steps leaving 
out technical details. Our purpose is to derive the qualitative structure of E(cr). To be explicit we choose harmonic 
interactions: 

Vl(x) = i ^ V nmX n X m , V nn = 0. (62) 

with coupling coefficients V nm = V(R n — R m )- Using V\(x) from Eq. 1)62(1 . Eq. takes the form: 

- x n + x z n - h = -e ^ VnnX m (63) 
m 

which can be solved by iterations: 

-x^ + {x^f -h = -eY J V nm x { ^ , x<® =x m (e = ) h)=x am (h) (64) 

I a 

etc. We remind the reader that x an (h) are the roots of Eq. J33J|. If \h — e^ m V nm x (Jrn (h)\ < h c for all cr, Eq. l(6^|) 
has three real roots: 

x^ (e, cr) = x„ n (h - e ^ V nm x am (h)) (65) 

i n 

which will be expanded with respect to e: 

x<p(e,tr) = x an (h)-ex' an (h)J2v nm Xa m (h)+ (66) 

m 

+ 7j e2x <r„( h ) ^2 V mnV nm >Xa m (h)x am ,(h). 



d L x„(h) 



can be represented as follows: 



%^=I>^K (67) 

i=0 

where as^ is easily be expressed by the 1-th derivative of x a (h), er = +, 0, — . Using Eq. I|67|) we get from Eq. (|66|l 

2 2 

4 1] (e>o-) = 5Z^„W<"+e X! I «.*»i»( ,l ) <T >m + (68) 

2„=0 i„j m =0 

1 2 

with: 



2 



sc< n (h) =a?:°'(ft), $ n ,i n i m (ti) = —a C W VnmXj m (h) , etc. . (69) 

m 

The substitution of a4 (e, 0*) from Eq. (J55) into V(a?) with Vi(cc) from Eq. leads to: 

£«(«r) = V{x^{ ( T))=Y J [e ( v\t) + e { i\e)o n + e { l\e)ol]+ (70) 



4, 

"1^12 



"b 6 ^ ] [•^nin 2 Ti3 (7 ni ff H2 ff i3 "b ^nin 2 n 3 (J n 2 < J n 

+ CTni^ 2 Cr„ 3 +0- ni <T n2 CT^ 3 ) +C^ln 2 n 3 ( a ' 2 ii CT n2 Cr n3 + 

i 2 2 , 2 2 \ , n (l) 2 2 2 1 , I \ 

+ ^n 1 IJ n2^ n3 + (T n^n 2 (T n 3 ) + ^ninafia a n t a n 2 a n 3 l + V ") 
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where (• • •) are four-, five- etc. body interactions which are of order e 3 ,e 4 etc. The coefficients in capital letters 
in Eq. (|70|l can be expressed by Xi n (h), x n ,i n i m (h), etc., h and V nm . E^(cr) is a kind of generalized Blume- 
Emery-Griffiths model. Using higher order iterates the corresponding energy E^'(cr) will be similar to Eq. H7()|l 
with "renormalized" coefficients ef\ A^ n2 , etc. On the other hand it is obvious that any function /(er) can be 
represented by a form as given on the r.h.s. of Eq. i|7U|) . Therefore we choose for E(cr) the r.h.s. of Eq. I|7l)|l without 
the superscripts. 

We begin with the calculation of n s (v). Taking X*(v, n s (v)) = (cf. Eq. (|ltJ|> ^ into account the saddle point 
equations (|TT))l reduce to: 

' ^E^e-^' 1 *^ , (71) 



N ■ Z(0,/j,*;N) 



1 - > ,-i»*N s {<r) 



N- Z(0,fi*;N) 
H*(n s ) is easily obtained from the second equation of (|71|l : 

^(n s )=iln-^2- (72) 
1 - n s 

and the first one can be written as follows: 

C = ^{H[ e o(e) + e i(e)K>o + e 2 (e)<^)o]+ (73) 

n 

+ 6 E [^"in2( cr ri 1 )o(c r »i 2 )o + S„ 1 „ 2 (((T, 2 11 )o(c r ri 2 )o + ("ni )o (<^ 2 )o) + 

+ C nin2 {^ 1 )o(o'l 2 )o} H } 

with 

(/(<0>o = E /(^e-^ 1 ^/ £ e-*"^ 1 -"') . (74) 

CT=+,0,- CT=+,0,- 

Taking /z* from Eq. (|72[) into account one easily finds: 

(tr„)o = , (o£) = 1 - n, . (75) 

This result is obvious, since for N s — n s N fixed the maximum weight of expression l)45|l is obtained for N + = N_ = 
N(l — n s )/2 which immediately implies Eq. i|75)l. Introducing the result from Eq. (|75|) into Eq. (|731 yields: 

v(n s ) = e (e) + e 2 (e)(l - n a ) + e<5 2 (0)(l - n s f + e 2 C 3 (0)(l - n s ) 3 + • • ■ (76) 

with: 

C-a(O) = £ Cta , <%(0) = E ^On,*,,-" (77) 

n 2 #0 0#n 2 #n 3 #0 

where the lattice translational invariance has been taken into account. For e small enough, v(n s ) is monotonous in 
n s and we can solve Eq. I|7fi|) for n s (v): 

r _v 1 , ^-(e (e) + e 2 (e)) <7 2 (0) - (e (e) + e 2 (e)) 2 2 , 

n s (w) = l + — e— — ( — ) +0e). (78) 

e 2 (e) e 2 (e) e 2 (e) 

Note that the bare one-particle quantities ej (cf. Eqs. H30() . are "renormalized" to ej(e) which depend on 

the coupling constants V^ m . Putting in Eq. (|78^l e = 0, one recovers (by use of Eq. lj2U), the result J2J). The 
interaction between the particles has two effects. First, it "renormalizes" the one-particle coefficients ei,i = 0,1,2 
such that ei(e) ^ 0, in general. This corresponds to a non-symmetric, effective on-site potential. Second, n s (v) 
becomes nonlinear. Its curvature depends on the sign of the effective coupling constant C 2 (0). 

The calculation of n s (v) is much more involved, because one has to perform averages with respect to exp[— iXE(<r)], 
which contains the interactions. This is in contrast to the averaging with exp[— i^iN s (cr)] (cf. Eq. (|7 1|) ^ which fac- 
torizes. Therefore an analytically exact determination of n s (v) for all v is not possible. But, in the following we 
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will demonstrate that the "thermodynamic" formalism can be applied to get fi s (v) for v/v gs - 1 < 1. In that limit 
we can perform a cumulant expansion. The simplest case is a "ferro-elastic" ground state, i.e. a^ s = + , as for the 
non-interacting particles. For v close to v gs = E(+, . . . ,+)/N only such stationary points exist for which a low 
concentration of a' n s deviate from +. Therefore we introduce "defect" variables: 

r n = 1 - cr n . (79) 

The ground state belongs to r„ = 0. r n = 1 or 2 indicates a "defect". Replacing in Eq. (|70|l a n by r„ yields: 

E(t) = E(+,...,+)+E (t)+E 1 (t) (80) 

with 

E (r) = ^[ei(e)r„+^( e )^] (81) 

n 

E i( T ) = ^2 [A nin2 (e)T ni T n2 + B nin2 (e){T 2 ni T n2 + r ni r^ 2 ) + Cn^i^T^r^} + 

+ ... 

where A nin21 B ni7l2 , etc. contain different orders in e and vanish for e = 0. They can easily be expressed by 
B nin2 , etc. For ei(e),e(e) one obtains: 

ei(e) + e 2 (e) = e (e) - v gs (e) = v (e) - v gs (e) (82) 
2(ei(e) + 2e 2 (e)) = e (e) - e x {e) + e 2 (e) - u ss (e) = u_(e) - v gs (e) 

where Vq(c) and u_ (e) are obtained from Eq. (|31|l by replacing by ej(e). Using a cumulant expansion we get for 
the "free energy" : 

/(A,/i) = % S A + / (A,m) + A^( j Bi(t)) + (83) 
+zA 2 ^[((i? 1 (r)) 2 >o-((i? 1 (r)» 2 ] + ... 
The quantities with subscript are obtained with the unperturbed "canonical ensemble" : 

po ( T ) = 1 -i(«>(r)+^.(r)) ( 84) 

where N s (t) is obtained from Eq. IMl'l) by use of Eq. I|79l) . Then the saddle point equations ll(j|) are as follows: 

v-v gs = ^ + 1^E 1 (t)) + X^(E 1 (t)) )+- (85) 

n * = ^7 + aV lM>0 + ^ ' 

Here, we have only given terms up to the first cumulant. /o(A,/i) can easily be calculated from Eq. I|84|) . and J£ 
are linear in (r„) and (r 2 ) where: 

A + 2B , 2v A + 4B 
{Tn)o = 1 + A + B ' (T " )o = 1 + A + B (86) 

with: 

4 = exp[-i(A(ei + e 2 ) + /i), ], 5 = exp[-i2A(ei + 2e 2 )]. (87) 

The correction terms in Eq. 1(8 5|1 . i.e. the first, second, etc. cumulants are quadratic, cubic, etc. in (r n )o and (t 2 )o- 
The l.h.s. of Eq. 1(85(1 becomes arbitrary small since v — v gs — » implies rt s — >■ 0. Then it follows that (r n )o and (t 2 )o 
in leading order are linear in v — v gs and n s and that the cumulant terms in Eq. i(85(l are of higher order. After having 
performed in Eq. 185(1 the differentiations, we are allowed to set fi(v,n s ) = 0, since we are calculating n s as function 
of v. Writing A and B as the result for e = plus a correction: 

A ~ ft , X \ R (V-Vgs) - (gl +e 2 )» 3 

= ns(1 + ^ ' B = 2(e 1 + 2e 2 ) (1 + fe) (88) 



13 



one finds from Eq. (|85() : 



S A ~ (an s + b(v - w ffs )) lnn s , — (cn s + d(v - v gs )) lnn s 



(89) 




with 



a; = 



(91) 



«_(e) -u+(e) 



and 



8{e) 



^o(e) -^ gs (e) 
i>_(e) - u 5S (e) 



(92) 



where we used Eq. 1)82(1. Comparison of the result l|9l)|l with the corresponding one for noninteracting particles 
(Eq. H43fl ) shows that the interactions do not change the power law dependence. But, they lead to modified next-to- 
leading order corrections, containing logarithmic dependence on v — v gs . 



The purpose of the present paper has been twofold. First of all, we wanted to establish a relationship between 
the stationary points which are solutions of a set of coupled, nonlinear equations, and symbolic sequences er, where 
o~ n takes values from an alphabet A. We have proven for the class of 4 -models that such a unique relation exists, 
provided that the coupling parameter e is below a critical value e c . This proof is based on Aubry's anti-continuum 
limiU2il& and the application of the implicit function theorem, similar to the proof of existence of breather by MacKay 
and Aubry 19 . Assuming that such a one-to-one correspondence exists in general we have demonstrated in the second 
section that one can use a "thermodynamic" formalism for the calculation of statistical properties of the stationary 
points. Consequently, this description allows to investigate topological quantities of the so-called "energy landscape" 
of a potential energy V(xi, . . . ,x n ). Particularly, it yields an explanation why a topological singularity, e.g. for the 
Euler characteristic of the manifold Ai n (v) (cf. Eq. can occur. We stress that our "thermodynamic" formalism is 
not the same as that recently used to calculate the full canonical partition function^. There the configuration space 
has been divided into basins of the stationary points. Taking also the vibrational degrees of freedom into account 
these authors were able to calculate the canonical partition function, under a couple of assumptions. 

The relationship between stationary points and sequences er already proves that there arc in total = cxp(N In nX) 
stationary points, ua is the number of "letters" of the "alphabet" A. For rather special models this has already 
been shown to be truei^2^24. Then it is easy to calculate the saddle index distribution function Ppf(n s ) because the 
number of negative eigenvalues of the Hessian can be related to the number of cr n 's taking certain values. For the 
</> 4 -model this value is 0. Simple combinatorics leads for P^(n s ) to a Gaussian of width proportional to N' 1 ^ 2 and 
a maximum at n™ ax = 1/3. It is interesting that this value obtained for a 4 -lattice model coincides with the result 
found for Lennard-Jones clusters for 4 < N < 9 3 . Is this just an accident? The answer is not clear. It may be true 
that the potential energy landscape of a liquid can be decomposed into basic "units" which are double-well-like. Of 
course, the smallest "unit" one can choose are two adjacent local minima. Because both must be connected by a 
barrier one arrives at a double- well potential with extrema being labelled by +, 0, — . The main open question is: Can 
one really construct the full energy landscape by connecting such double well potentials and accounting correctly for 
its connectivity? 

For the 4 -model we have found the energy dependence of the saddle index. It has become clear that the functions 
n s (v) and n s (v) are not identical. Whereas both functions are equal one at the largest possible value for v and v 
■Vmax — v m ax = iV _1 .E(0, . . . , 0), they have different behavior below N~ 1 E(0, . . . , 0). n s (v) vanishes at v* above the 
ground state energy v gs and n s (v) becomes zero at v = v gs , only. We think that this behavior is true even for liquid 
systems. The "thermodynamic" formalism has made it possible to get quantitative results for both functions, at least 
for small enough coupling parameter e. n s (v) is linear in v for e = and becomes nonlinear for e > 0. n s (v) is already 
nonlinear in v in case of e = 0. For v close to v gs it exhibits a power law. We have shown that this power law behavior 
is not changed for < e < e c , assuming the ground state to be ferroelastic. For an antiferroelastic ground state one 
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arrives at the same conclusion. However, whether the power law exists for periodic ground states with period larger 
than two or even for quasiperiodic ones is not yet clear. In addition, its importance for physically relevant quantities 
is not obvious. 

The fact that functions n s and n s are different, seems to be in variance with the numerical results of Ref^. These 
two functions were found to be the same, within statistical errors. This implies that also n s (v) vanishes at an energy 
above v gs . The liquid in Rcf.« was equilibrated at a temperature, e.g. T — 2 and 0.5 (in Lennard-Jones units). T = 0.5 
is already close to the mode coupling temperature T c where the dynamics is already rather slow. For n s = there 
is a huge variety of stationary points with energies v gs < v < N~ 1 E(cr) with a n = + or — . Their maximum weight 
occurs when N + — N- = N/2 (N a = number of er n 's in er which are equal to a). The stationary points with n s = 
which are close to v gs have an exponentially smaller weight and can only be found at T much smaller than 0.5 where 
the configurations look like a crystal with a low concentration of defects. It seems to us unlikely that the simulation 
done in Refi^ which has reached T = 0.5 from the liquid phase has really been in the range of such defected crystalline 
configurations. Our result that n s (v) vanishes at v gs , only, which implies that n s as function of temperature vanishes 
at T = 0, only, is consistent with recent results (— and second paper of&). 

Let us come back to a liquid system. How could one apply a similar strategy as used in the present paper? The 
answer is as follows: Divide the sample with mN particles into m-boxes with N particles where 1 <C m <C N. Let us 
switch off the inter-box interactions. Then the number of stationary points and the saddle index properties can be 
related to the corresponding quantities of a single box. This starting point corresponds to Aubry's anti-continuum 
limit, although a single box represents already a nontrivial problem. Then one can use again the implicit function 
theorem to prove that the topological features are unchanged under turning on the inter-box interactions, provided 
their strength is below a critical value. However, increasing this strength such that inter- and intra-box- interactions 
are the same may exceed the critical strength. As far as we know this type of reasoning has been used first by 
Stillinger (see2^ and references therein). Recently it has been used again. 22,26 . Based on independent boxes the 
authors of Ref^& have derived a relationship between a, a and 7 which yield the number of stationary points with 
saddle index n s = N s /N of a single box: 

M n (n s = 0) = expaN (93) 
M n (n s = l/N) = aN expaN 

and the total number of stationary points of a box: 

M^ t = exp 1 N. (94) 

This relation is: 

7 = a + aln2. (95) 
Now let us apply this reasoning to our -model. There it is: 

M N (n s = 0) = 2 N = exp(Aln2) (96) 
M N (n s = l/N) = N ■ 2 1 *- 1 = iiVcxp(iVln2) 
Af*f = 3^ = exp(iVln3) , 
i.e. we find: a = ln2, a = 4, and 7 = In 3. Substituting a and a into Eq. (|95|l we get for its r.h.s.: 

Ir.h.s. = ^ln2~ 1.039 (97) 

which is close but not identical to 7 = In 3 ~ 1.098. Whether this small discrepancy is a hint that correlations between 
boxes can be neglected or not, is not obvious. 

To conclude, we have demonstrated the usefulness of a one-to-one correspondence between stationary points and 
symbolic sequences. 

Acknowledgement: I am very grateful for all the discussions I have had with Serge Aubry. The present paper 
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